Accepted for publication in the Astrophysical Journal 

Preprint typeset using I^T^X style cmulatcapj v. 03/07/07 



O 

o 

(N 



00 
(N 

o 

u 

Oh- 
i < 

o ■ 



> 
oo 
m 

O 

&\ 
O 



13 



STRONG LENSING BY SUBHALOS IN THE DWARF-GALAXY MASS RANGE II: DETECTION 

PROBABILITIES 

T. Riehm 1 , E. ZACKRISSON 2 ' 1 ' 3 , E. Mortsell 4 & K. Wiik 2 
Accepted for publication in the Astrophysical Journal 

ABSTRACT 

The dark halo substructures predicted by current cold dark matter simulations may in principle be 
detectable through strong-lensing image splitting of quasars on small angular scales (0.01 arcseconds 
or below). Here, we estimate the overall probabilities for lensing by substructures in a host halo closely 
aligned to the line of sight to a background quasar. Under the assumption that the quasar can be 
approximated as a point source, the optical depth for strong gravitational lensing by subhalos typically 
turns out to be very small (r < 0.01), contrary to previous claims. We therefore conclude that it 
is currently not feasible to use this strategy to put the simulation predictions for the dark matter 
subhalo population to the test. However, if one assumes the source to be spatially extended, as is 
the case for a quasar observed at radio wavelengths, there is a reasonable probability for witnessing 
substructure lensing effects even at rather large projected distances from the host galaxy, provided 
that the angular resolution is sufficient. While multiply-imaged, radio-loud quasars would be the best 
targets for unambiguously detecting dark matter subhalos, even singly-imaged radio quasars might 
be useful for setting upper limits on the abundance and central surface mass density of subhalos. 
Subject headings: gravitational lensing — galaxies: halos — quasars: general — dark matter 



1. INTRODUCTION 

In the standard cold dark matter (CDM) scenario, as 
well as in several slight modifications thereof, dark mat- 
ter halos are assembled hierarchically from smaller sub- 
units. At the time of merging, some of these subunits 
are disrupted and dispersed into the smooth dark mat- 
ter component of the halo, whereas others temporar- 
ily survive in the form of subhalos. A long-standing 
problem with this picture is that the number of sub- 
halos predicted by CDM simulations is orders of mag- 
nitudes higher than the known number of s atellite galax- 
ies in the vicinity of the Milky Way (e.g. iKlvpin et al.l 
Il999t iMoore et alj|199fl iDiemand et a 1.1 k200"l J). There 
are several possible ways out of this dilemma: either the 
CDM scenario is incorrect, the simulation predictions are 
wrong, or the majority of these subhalos must somehow 
have evaded detection. The latter alternative is quite 
viable, provided that the baryonic content of these sub- 
halos have been either lost or prevented to form stars 
dBullock et al.ll2000tlB"enson et al.ll2002l:ISomervillejl2002l: 
Kravtsov et all l2004l ~ lMoore et all 120061 : IStrigari et al I 
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M© range would be gravitational millilens- 
image splitting a t a characteristic separa- 
tion of milliarcseconds (e.g. IWambsganss fc Paczvnskil 



1992t iBarvshev fc Ezoval Il997t iMetcalf fc Madaul l2001fc 
Yonehara et al. 1 12003ft . Based on a null detection of 



2007T ) . Recent studies have shown that when including 



these physical processes together with observational bi- 
ases most of the predicted Milky Way sat ellites would lie 
outside the reach of curren t surveys fe.g iTollerud et alJ 
I2008t IKoposov et all 120091: IMaccio et al.ll2009fh If such 
"dark galaxies" do indeed exist, gravitational lensing 
may offer one of the most promising ways to detect them. 
One tell-tale signature of dark matter subhalos in the 

Electronic address: teresa@astro.su.se 

1 Department of Astronomy, Stockholm University, AlbaNova 
University Center, SE-106 91 Stockholm, Sweden 

2 Tuorla Observatory, Department of Physics and Astronomy, 
University of Turku, Vaisalantie 20, FI-215 00 Piikkio, Finland 

3 Department of Astronomy and Space Physics, Box 515, SE-751 
20 Uppsala, Sweden 

4 Department of Physics, Stockholm University, AlbaNova Uni- 
versity Center, SE-106 91 Stockholm, Sweden 



millilensing in a sample of 3 00 qu asars observed with 
the VLBI, IWilkinson et al.l (|2001h demonstrated that 
the vast majority of quasars do not show any signs of 
millilensing, and were able to impose an upper limit 
of SI < 0.01 on the cosmological density of point- 
mass dark matter objects in this mass range. Un- 
fortunately, this limit is insufficient to set any use- 
ful constraints on subhalos predicted by CDM, since 
their lensing properties are very diff erent from those of 
point-masses (jZackrisson et al.l(2008l hereafter paper I). 
Adopting more realistic subhalo density profiles would 
significantly raise the above limit. To put the CDM 
subhalo predictions to the test, it has instead been sug- 
gested that one should target quasars which are al- 
ready known to be gravitationally lensed on arcsecond 
scales, as one can then be sure that there is a massive 
halo well-aligned with the line of sight, which substan- 
tially increases the pr obability for subhalo millilensing 
(| Yonehara et al.l l2003f h Indeed, the magnification as- 
sociated with millilensing has long been suspected to 
be the cause of the flux ratio anomalies seen in such 
syste ms (e.g. iMao fe Schnei der 1998; Ko chanek fc Dalai 
2004). However, current CDM simulations appear to 
predict too few subhalos to explain the flux ratios of 
many of these systems (e.g. IMao et al.l 2004 iMetcalJ 
l2005at IMaccio fc Mirandal 120061 : Ku et al.l 12009ft . Thus, 
other mechanisms like lensing by low-m ass field halos 
(jMetcalJ I2005ri iMiranda fc Macciol l2007ft. stellar mi- 
crolensing in the lens galaxy (jSchechter fc Wambsganssl 
2002) and absorption or sc intillation in the interstellar 
medium (jMittal et al.ll2007ft may also be at work. Sub- 
halo millilensing has also been advocated as a n expla- 
nation for strange bending angles of radio jets (|MetcalJ 
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2002) and image positions whic h smooth halo models 
seem unable to a ccount for (e.g. iKoopmans et alJl2002t 
iBiggs et al.M2004t iMore et al.ll2008[) . ~ 

In paper I, we estimated the image separations for the 
subhalo density profiles favoured by recent N-body sim- 
ulations, and compared these to the angular resolution 
of both existing and upcoming observational facilities. 
In this second paper in the series, we assess the proba- 
bility for subhalo millilensing of a quasar as a function 
of projected radius from the center of a foreground lens 
galaxy. The paper is structured as follows. In f[2] the 
properties of host halos are discussed. The subhalo pop- 
ulation is described in Sj3] In Sj4j we compute the optical 
depth for subhalo lensing of a point-source as a function 
of projected radius from the host halo. The expected 
number of intervening subhalos located in a galaxy-sized 
host halo along the line of sight to an extended source 
is estimated in Sj5] Finally, in ^ we discuss several ef- 
fects that could affect our predictions and present our 
conclusions. 

Throughout the paper, we assume a ACDM cosmology 
with n A = 0.762, n M = 0.238 and h = 0.73 (H = 100 h 
km s^ 1 Mpc " 1 ) in concordance with the WMAP 3-year 
data release (Spc rgel et al.ll200"7l ). 



2. HOST HALO PROPERTIES 

In this study, we focus on host halo lenses in the typical 
galaxy m ass range which span s 10 10 M Q < M < 2 x 10 13 
(e.g. ILTfe Ostriked 12001 . The upper limit is set 
by the fact that more massive dark matter halos cor- 
respond to groups or clusters of galaxies. Although it 
might be interesting to investigate the dark substruc- 
ture population in these systems, they are presumably 
fundamentally different from galaxy-sized halos. This is 
partly due to a difference in the distribution of the lu- 
minous component. Furthermore, it is only for a critical 
mass < 10 13 M© that cooling of the baryonic halo com- 
ponent can lead to a concentration of the baryons to the 
inner parts of the mass profile, resulting i n a steep in- 
ner density profile (jBlumenthal et aLlll986[ ). Therefore, 
halos at higher masses are expected to exhibit shallow 
inner density p rofiles similar to the profile proposed by 
iNavarro et alJ (|1997l hereafter NFW), w hich has been 
confi r med by several lensing s tudies (e.g iLi k, Ostrikerl 
l2002t iMandelbaum et al.ll2006f ). Most halos with masses 
below 10 10 M© fall in the dwarf galaxy mass range, which 
corresponds to the subhalo masses we are interested in. 
A host halo in this mass range would comprise subha- 
los with even smaller masses (typically lighter than 1% 
of the host mass, see below) and therefore not cover the 
entire relevant subhalo mass range. 

Unless stated otherwise, we assume the dark subhalos 
to reside within a Milky- Way like host halo of mass M = 
1.8 x 10 12 Mq at a redshift zi — 0.5, while the sources 
are placed at redshift z s = 2. 

The singular isothermal sphere (SIS), with density 
profile p oc r~ 2 has proved to be a successful model 
for luminous galaxy-mas s lenses (e.g. iRusin et alJ 120031 : 
iTreu & KoopmansH2004l) . This is believed to be due to 
the substantial contribution from baryons to the inner 
regions of these objects, as N-body simulations based on 
the CDM scenario predict the dark matter halo to be sub- 
stantially less centrally concentrated (NFW). The mass 



of an SIS inside a three-dimensional radius R is 

2o*R 



M(R) 



G 



(1) 



with the one-dimensional velocity dispersion a v . The 
projected surface mass density of an SIS is then given by 

„2 



2GC 
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where £ is the projected distance from the center of the 
lens. 

The convergence k and the shear 7 for an SIS are iden- 
tical and simply prop ortional to its surface density (e.g. 

umm-- 



ISchneider et alJ 



k = 7 = 



with the critical surface density defined as 

c 2 D x 



u crit 
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where D s , D\ and D\ s represent the angular diameter dis- 
tance from the observer to the source, from the observer 
to the lens, and from the lens to the source, respectively. 

The lensing magnification from the host halo /1 is then 
given by 



/i=[(l- K ) 2 - 7 2 



(1-2/i)- 1 . 



(5) 



The host halo is truncated at the virial radius -R200 
defined as the radius at which the mean enclosed density 
equals 200 times the mean mass density of the Universe 
at redshift z, 
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3. SUBHALO POPULATIONS 



We are interested in dark subhalos with masses m within 
the range 10 6 < m/M© < 10 10 . These limits are partly 
due to the current resolution of the numerical simulations 
used to predict the existence of substructure. Further- 
more, as has been shown in paper I, the lower mass limit 
roughly corresponds to a minimum halo mass for which 
the expected image separation could be resolved with ei- 
ther current or upcoming observational facilities. The 
upper mass limit roughly corresponds to the mass below 
which the discrepancy between the number densities of 
luminous galaxie s and dark matter halos starts to be- 
come severe (e.g. IVerde et alj|2002t Ivan den Bosch et"aL1 
12001 . 

We assum e a subhalo popu lation following the model 
proposed by iGao et ahl (|2004f ) . The subhalo abundance 
per unit halo mass (ignoring the high mass cut-off, m > 
0.01 M) can be approximated by 



^ = 10" 
dm 



3.2 



-1.9 



KM? 



(7) 
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iDiemand et al.l (|2007af ) confirmed this finding and ex- 
tended it to lower subhalo masses in a recent high reso- 
lution simulation of CDM substructure in a Milky Way- 
sized halo. A fraction f su b{M, m m i ni m max ) of the mass 
M of a host halo is in the form of subhalos with min- 
imum and maximum masses m m i n and m maxi respec- 
tively. Here we assume subhalo masses m in the range 
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4 x 10 6 < m/Mp, < 10 10 , corres ponding to the interval 
probed bv lDiemand et all ()2007al ). The total mass frac- 
tion of a host halo with M — 1.8 x 10 12 M Q in subhalos 
is then about 5% for this subhalo mass range. 

Within this simulation, it has also been shown that 
the subhalo nu mber density profile can be fitted by the 



tne subhalo nu mber density profile < 
following form (|Madau et al.ll2008f) : 



n(< x) 
N 



12x s 



1 + 11a; 2 : 



(8) 



where x is the distance to the host center in units of i?2oo, 
n(< x) is the number of subhalos within x and N is the 
total number of subhalos inside i?200- 

4. POINT-SOURCES 

In the optical wavelength region, a quasar can for the 
purpose of this paper be approximated by a point source. 
Under the assumption that the lenses do not overlap 
along the line of sight, the optical depth r represents 
the fraction of a given patch of the sky that is covered 
by regions in which a point source will be lensed. 



T"(£»«,7) = -g 



<7iens(m, K,j) dl ]^ dm (9) 
dm 



with the area of a patch on the sky S and n(£) the number 
of subhalos projected on S at the projected distance £ 
from the center of the host halo which has been derived 
numerically using equation ([8]). Here, k and 7 are the 
external convergence and shear induced by the host halo, 
respectively, and ai ens denotes the cross section for a 
single subhalo lens within the potential of the host halo. 
In the limit of small r, the optical depth can therefore 
d irectly b e used as an estimate of the lensing probability. 

iKeetonl (|2003h derived an analytic expression for the 
lensing cross section of an SIS subhalo within a host halo 
potential, 



<7i ens (m, 5, re, 7) = fxA(m,S,K,j), 



(10) 



where A corresponds to the area in the source plane 
where the total magnification perturbation due to the 
substructure is stronger than 6. In the limit \5\ — > 00, this 
analysis gives the area inside the caustic where we expect 
to obtain multiple images from substructure millilensing. 
A discussion on the validity of the SIS assumption for 
substructures can be found in [JSJ The area A can be 
expressed through the Einstein radius of the subhalo, 6, 
which is defined as 



4"7T 



V c J D s 



(11) 



with the conventional conversion between mass m and 
the subhalo velocity dispersion a v : 



Gm 
2r 20 o ' 



For \S\ — > 00, A is given by 
A(m,K,j) 



-irb 2 j 2 \x 



(12) 



(13) 



We can use this formalism to compute the optical 
depth for subhalo lensing as a function of projected dis- 
tance £ from the center of the host halo. At the Einstein 



radius of the host halo, the magnification diverges to in- 
finity which causes the optical depth at this radius to 
diverge as well. However, this is not physical since even 
for point sources the maximum magnifi cation is limited 
due t o the effects of wave optics (e.g. ISchneider et al.l 
1992). Furthermore, the area in the source plane for 
which a background source would experience very high 
magnifications is rather small. Of the approximately 1 in 
500 quasars which are strongly lensed, only around 1 in 
(/i mM — l) 2 would experience magnifications exceeding 
Mmax- This illustrates the unlikeliness of a galaxy lens 
producing very high magnifications on a random back- 
ground quasar (circa one per 5 million for /i max = 100). 
Of course, any given sample of multiply imaged quasars 
will be biased toward the high magnification tail since 
systems with low magnification will fall short of the mag- 
nitude limit of the survey. The size of this magnification 
bias will depend on the quasar luminosity function (QLF) 
and the magnitude limit of the survey. In recent years, 
there has been a great effort to constrain the QLF us- 
ing various large surv eys producing thousands of QSOs 
(e.g.lBovle et al.ll2000t[Croom et alj2004HRichards et all 
I2006al) . However, it has been proven rather difficult to 
trace the QLF simultaneously to both faint magnitudes 
and hig h redshift. Here w e use the QLF for z ~ 3.2 
given in ISiana et "all (|2008l ) covering QSOs at faint-end 
magnitudes and convolve it with the lensing cross sec- 
tion for a singular isothermal ellipsoid to estimate the 
effects of magnification bias. For the QLF magnitude 
limit r' < 22, we compute that approximately 0.4% of 
all observed multiply-imaged QSOs will experience mag- 
nifications // > 100 (compared to ~ 0.01 % from our sim- 
ple estimate not taking magnification bias into account). 
Lowering the magnitude limit to r' < 20 roughly doubles 
the fraction of multiply lensed quasars with fi > 100, 
thereby demonstrating the importance of the magnitude 
limit. Although we expect z ~ 3.2 to be a typical value 
for an observed source redshift, we are also interested 
in higher rcdshifts where the effects of magnifaction bias 
will b ecome even more pronounced. In IWyithe fc Loebl 
(2002), it was shown that for the Sloan Digital Sky Sur- 
vey (SDSS) with magnitude limit i* < 20 (z* < 20.2), 
the fraction of QSOs which are expected to be multiply 
imaged at z ~ 4.3 (6.0) might be as high as ~ 13 % (30 %) 
of which ~ 5 % (10 %) should experience magnifications 
jU > 100. However, recent results from the SDSS suggest 
that the QLF bright-end slo pe at z > 3 is getting shal- 
lower toward higher redshift (jRichards et al.ll2~006al ) and 
high-resolution HST observations of 157 SDSS QSOs at 
4.0 < z < 5.4 resulted in the nondetection of strong lens- 
ing in these systems (jRichards et al.l l2~006bf ) . Thus, we 
put an upper limit /i max = 100 for a realistic expected 
maximum magn i ficatio n from the host halo and use the 
IWyithe fc Loebl (|2002l ) formalism as a conservative esti- 
mate of the effect of magnification bias at high redshift. 
Currently, the highest redshift quasar with multiple im- 
ages in the SDSS is at z s = 3.626 w ith its lensing galaxy 
at zi = 0.4 - 0.6 (jlnada et al.ll2008[ ) while the most dis- 
tant multiply-imaged quasar kno wn today is at z ? = 4. 5 
with a lensing galaxy at zi = 0.6 (McMah on et al.| [l992). 

As can be seen in figure HJ the optical depth for 
Mmax = 100 does not exceed a value of 0.005 at any ra- 
dius. This is contrary to prior claims where the optical 
depth for subhalo lensing had been estimated to be orders 
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Fig. 1. — Optical depth for subhalo millilensing of a point source 
at z s = 2 as a function of projected radius from the host halo 
center. Here we assume a host halo at z; = 0.5 with mass M = 
1.8 X 1O 12 M and subhalos in the mass range 4 X 10 6 - 10 10 M Q . 
The peak has been cut with respect to a maximum magnification 
factor of 100 at the Einstein radius of the host halo. 




Fig. 2. — Optical depth for subhalo millilensing corresponding 
to a maximum magnification factor of 100 at the Einstein radius 
of the host halo for a point source as a function of lens redshift, 
zi, and source redshift, z s . We assume a host halo with mass 
M = 1.8 X 10 12 Mq and subhalos in the mass range 4 X 10 6 - 
1O 1O M . 

of magnitude higher. lYonehara et alj (|2003f ) evaluated 
the quasar millilensing optical depth from SIS subhalos 
with a slightly tighter mass interval (10 7 - 10 10 M ) to 
be approximately 0.1. If we reconstruct their scenario, 
we get optical depths below 10~ 3 . This discrepancy can 
be traced to an incorrect subhalo mass function adopted 
bv lYonehara et alj (|2003t) . 

Figure [2] shows the dependence of the optical depth r 
on the lens and source redshifts, zi and z s respectively. 
Here we plot the optical depth at the Einstein radius of 
the host halo lens with a maximum magnification (i max 
— 100. These optical depths arc therefore upper lim- 
its which will only be valid for quasars favorably aligned 
with the host halo. Since the optical depth for subhalo 
lensing increases with both lens and source redshift, ob- 
serving high-z objects will increase the probability for 
detecting such an effect considerably. Still, it would re- 
quire an immense amount of fine-tuning concerning the 
lens-source alignment in both projected radius and red- 
shift to reach optical depths exceeding 0.025. 

In order to investigate how these results might change 
taking magnification bias into account, we also esti- 
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Fig. 3. — Optical depth for subhalo millilensing of a point source 
at z s = 2 as a function of projected radius from the center of 
a host halo at z\ = 0.5 with mass M = 10 1:l Mq (dotted line), 
1O 12 M0 (dashed line) and 1O 13 M0 (solid line), respectively. For 
M = 1O 11 M0, rrimax has been set to 10 9 Mq. The peak has been 
cut with respect to a maximum magnification factor of 100 at the 
Einstein radius of the host halo. 

mate the expectation value for the optical depth < r > 
by relaxing the concept of the maximum magnification 
[i-max and instead using the magnification distribution for 
mult iply- imaged quasa rs in the SDSS given in figure 6 in 
IWvithe fe Loebl (|2002ft . The trend for higher lens and 
source redshifts to give rise to higher expectation values 
for the optical depth < r > is amplified by the effect of 
magnification bias. While the expectation value for the 
optical depth for a system with z\ = 0.5 and z s « 2 is 

< t >^ 10 -4 , the corresponding values for z s 4 and 6 
are < r >^ 10~ 3 and 4 x 10~ 3 , respectively. Raising the 
lens redshift for the latter two scenarios to zi — 2, results 
in an additional increase of the expected optical depth 

< t > by less than one order of magnitude, only reach- 
ing values below 0.05. However, as ment ioned above, the 
effect of magnification bias described in IWvithe fc Loebl 
(2002) has been shown to be conservative and there are 
no systems with such redshift combinations known to- 
day. Thus, only a dedicated high-z survey, collecting a 
substantial number (> 100) of such objects, might chal- 
lenge our conclusions concerning the bleak prospects for 
subhalo detection through quasar image splitting in the 
optical. 

In figure [31 we explore the dependence of r on the 
host halo mass. For a point source at z s — 2 the opti- 
cal depth for subhalo lensing as a function of projected 
radius £ has been computed for a host halo at zi = 0.5 
with mass M = 10 11 M , 10 12 M and 10 13 M , respec- 
tively. For a host halo mass of 10 11 M we limit m ma x to 
10 9 M since the subhalo mass function is not valid for 
m > 0.01A/. Subhalos above this mass limit would also 
resemble a galaxy-galaxy merger during accretion to its 
host halo and disturb the system, thereby invalidating 
previous assumptions. Larger host halos contain a larger 
fraction in substructure which results in an increased op- 
tical depth. However, this effect is not strong enough to 
substantially raise the probability for detecting subhalo 
lensing. Even for a host halo mass of 2 x 10 13 M which 
is our upper mass limit for a galaxy-sized host halo, the 
maximum optical depth lies below 0.01. 

In summary, we find that the optical depth for strong 
gravitational lensing of point sources by subhalos is 
lower than previously assumed, typically well below 0.01. 
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Fig. 4. — Average number of substructures projected on an 
extended source at z s = 2 as a function of projected radius for 
different source sizes. We assume a host halo at z; = 0.5 with a 
mass M = 1.8 X 1O 12 M0 and subhalos in the mass range 4 X 10 6 
- 10 10 Mq. We plot our results for a source with radius r s = 10 pc 
(dotted line), 100 pc (dash-dotted line), 1 kpc (dashed line) and 10 
kpc (solid line) , respectively. We assume a maximum magnification 
factor of 30 at the Einstein radius of the host halo. 

Hence, unless a substantial sample of multiply- imaged 
quasars would become available, a search for quasar 
image-splittings by subhalos is unlikely to result in any 
detections. 

5. EXTENDED SOURCES 

If one considers observing the subhalo lensing effects on a 
quasar at longer wavelength, the quasar can no longer be 
approximated by a point source. For extended sources, 
the lensing effects would appear as monopole-like or 
dipole-like distortion patterns in the surface bright- 
ness profile of the sour ce rather than multiple imaging. 
Ilnoue &: Chibal (|2005bh have argued that it should be 
possible to detect these lensing effects from dark matter 
substructures when observing extended sources resolved 
at scales smaller than the Einstein radii of the subha- 
los. Furthermore, it should even be possible to put con- 
straints on the internal density profiles of the lensing sub- 
halos. This technique may already become obser vation- 
ally feasible with ALMA 5 (|Inoue fc Chibal l2005al) or fu- 
ture space- VLBI missions like VSOP-2 B (jlnoue fc Chibal 
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l2005d ) and thus constitute a major step forward in the 
study of dark halo substructures. Here we estimate the 
probability of subhalo lensing in the radio regime by com- 
puting the average number of subhalos that lie within the 
region of the host h alo covering the source. 

In a recent study, iTorniainen et al.l (|2008| ) have shown 
that the typical source size for a quasar ranges from sev- 
eral 10 pc to a few kpc when its turnover frequency falls 
in the radio regime. 

As has been shown bv lPerrotta et all (|2002f ). the max- 
imum magnification by a galaxy-sized halo that can be 
achieved for extended sources with an effective radius of 
1-10 kpc at redshifts within z = 1 - 4 falls into the 
range 10 - 30. We therefore set a conservative upper 
limit of 30 for the maximum magnification from the host 
halo potential. 

In figure[H the expected number of subhalos projected 
on an extended source at z s = 2 as a function of pro- 
jected distance from the host halo lens center is shown 

5 www.alma.info 

6 www.vsop.isas.ac.jp/vsop2 
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Fig. 5. — Average number of substructures projected on an 
extended source with radius r g = 1 kpc at a projected radius of 
R200/2 (ranging from 22-7 arcseconds for z; = 0.5 - 2) as a 
function of lens redshift, z;, and source redshift, z s . We assume a 
host halo with mass M = 1.8 X 10 12 Mq and subhalos in the mass 
range 4 x 10 6 - 1O 1O M . 



for several source radii r s , ranging from 10 pc to 10 kpc. 
This can be compared to the virial radius -R200 ~ 260 
kpc for the host galaxy at zi — 0.5. It becomes clear 
that for sufficient source size (> 1 kpc) there is a good 
probability for the source image to be affected by subhalo 
lensing, not only close to the Einstein radius of the host 
halo but even at a rather large projected distance from 
the host halo lens center. For a source with r s = 1 kpc, 
one would expect at least one intervening subhalo per 10 
observed systems with a maximum projected distance 
of 10 arcseconds between the foreground galaxy and the 
source. For r s = 10 kpc, this number increases to ap- 
proximately 10 subhalos projected on the source out to a 
distance of 10 arcseconds from the host galaxy. This im- 
plies that for extended sources not only multiple-image 
systems but even the much more common singly-imaged 
quasars should be affected by subhalo lensing. 

Since quasars in the radio are likely to posses intrin- 
sic structure, one has to be careful not to confuse this 
with a potential lensing signal from subhalo lensing. As 
internal structures should be mapped on all images of 
a multiply-imaged source, distortions due to additional 
lensing of one image on small scales are distinguishable. 
However, even for singly-imaged quasar one can compare 
the observed amount of source structure to what would 
be expected from millilensing by subhalos low-mass field 
halos. Thus, it opens up the possibility to set upper lim- 
its on the abundance and central surface mass density of 
subhalos using singly-imaged quasars with a foreground 
galaxy projected within some tens of arcseconds of the 
source. 

Similar to the case of a point source, we also explore the 
dependence of these results on the lens and source red- 
shifts, zi and z s respectively. Since this scenario seems 
to be much more feasible than that of a point source, 
we are not only interested in upper limits but a realis- 
tic estimate. The probability for a quasar to be aligned 
close to the Einstein radius of a foreground galaxy, where 
its magnification will boost the expected number of in- 
tervening subhalo lenses, is low even when considering 
magnification bias. However, it is expected that there is 
a large number of galaxy-quasar pairs with a projected 
separation smaller than the virial radius of the host halo 
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Fig. 6. — Average number of substructures projected on an 
extended source with radius r s = 1 kpc at z s = 2 as a function of 
projected radius. Here we assume a host halo at z\ = 0.5 with mass 
M = 10 n M Q (dotted line), 1O 12 M (dashed line) and 1O 13 M 
(solid line), respectively. For M = 1O 11 M0, m max has been set to 
10 9 Mq. We assume a maximum magnification factor of 30 at the 
Einstein radius of the host halo. 

-R200 (typically up to a few tens of arcseconds). In figure 
the average number of intervening subhalos at i?20o/2 
for a source size of r s = 1 kpc is shown as a function 
of zi and z s . Also here, the probability for lensing by 
subhalos increases by up to one order of magnitude for 
high-z objects. 

Figure [6] shows the dependence on the host halo mass 
of the average number of subhalos intervening with an 
extended source. We show results for a host halo at 
zi = 0.5 with mass M = 10 11 M Q , 10 12 M Q and 10 13 
M respectively, and a source with radius r s = 1 kpc 
at z s = 2. As in the case for point masses, we limit the 
maximum subhalo mass m max to 10 9 M© for a host halo 
mass of 10 11 M Q . As expected, there is an increase with 
host mass in the average number of subhalos projected 
on the background source due to the higher fraction in 
substructure for high-mass host halos. At half the host 
halo virial radius -R200/2, this can boost the expected 
number of subhalos projected on the background source 
by one order of magnitude over the mass range of galaxy- 
sized host halos. In addition, high-mass host halos pos- 
sess larger virial radii -R200, allowing for larger projected 
distances between foreground galaxies and background 
sources. 

All estimates on the number of intervening substruc- 
tures, N, quoted above, however, strongly depend on the 
minimum subhalo mass, m m i n , set. Since the effects pro- 
duced by the smallest substructures might not be strong 
enough to be resolvable with any available or future in- 
strument, there might be an effective minimum subhalo 
mass above the limit considered here. Implications from 
this are discussed in ff5J 

6. DISCUSSION AND CONCLUSIONS 

We have shown that the optical depth r for subhalo lens- 
ing of point sources is lower than previously predicted. 
Even for favorable conditions the upper limit on r typ- 
ically lies well below 0.01. The highest optical depths 
are reached close to the Einstein radius of the host halo 
where the background source is expected to be multiply- 
imaged. This implies that for the best probability to 
observe any effects from subhalo lensing one would have 
to target multiply-imaged quasar systems at favorable 



redshifts. Today there are only around 100 of these sys- 
tems known 7 and essentially all measured lens redshifts 
lie below 1. Adopting the numbers for r as shown in 
figure [2j where a high host halo magnification of 100 is 
assumed, we expect that on average 2 images (and no 
more than 4 images at the 95% C.L.) out of 100 in this 
sample could show signs of dark subhalo lensing. Only 
for high-redshift multiply-imaged quasar systems with 
lens and source redshifts zi w 2 and z s > 4, respec- 
tively, the expectation value for the optical depth < r > 
reaches values close to 0.05 when including magnification 
bias. However, since no such systems are known today, 
an extensive high-redshift survey would be needed. 

Furthermore, one has to be careful when analysing 
millilensing signals, since they might not be attributed to 
halo substructure but low-mass field halos along the line 
of sight. Although the lensing efficiency for dark matter 
subst ructures peaks if they are associated with the host 
halo, iKeetonl (|2003f ) has shown that SIS substructures 
may be moved in redshift by several tenths and still have 
a significant lensing effect. 

It is also important to point out that our estimates are 
valid for substructures with SIS profiles. It has been 
shown that the Einstein radiu s of a lens strongly de- 
pends on its density pr ofile (e.g. lWright fc Brainerdl2000l : 
Zackrisso n et al.ll2008l ). Therefore, adopting a different 
density profile for these lenses will alter their lensing 
cross sections and thereby the expected optical depth. 
N-body simulations based on the CDM paradigm typi- 
cally predict dark matter halos to have inner density pro- 
files of the form p{r) oc r~ a with central density slopes 
awl (e.g. NFW), compared to the SIS with a = 2. 
The Einstein radius for an NFW type subhalo is sev- 
eral orders of magnitude lower than that for an SIS halo 
of the same mass. Assuming the subhalos density pro- 
files to follow that proposed by NFW would therefore 
lower the optical depth for point sources considerably. 
Eve n adopting a more favorable density profile with a = 
1.5 ([Moore et al.lll999l . hereafter M99) would lower the 
optic al depth computed substantially (| Zackrisson et ail 
2008) . Recent high-resolution simulations have indicated 
that substructures might have an inner slope slightly 
steeper than NF W with significant hal o-to-halo varia- 
tions with a ~ 1.2 (Die mand et al.ll2008f) whereas others 
favour even shallow er inner slopes (jNavarro et~aT][200l 
ISpringel et al.1 [20081 ). We conclude that it is currently 
not feasible to use this technique to search for strong 
lensing signatures of quasars in the optical. 

If one instead targets the radio wavelength regime 
where quasars appear as extended sources, there is a 
high probability for subhalo lensing of quasars of suf- 
ficient size. For source sizes r s > 1 kpc, this is valid 
even at rather large projected distance of the source to 
the host halo center. This allows for a different search 
strategy than those previously proposed. Instead of only 
targeting multiply-imaged quasar systems, even quasar- 
galaxy pairs with a separation of several tens of arcsec- 
onds should show effects of strong lensing by substruc- 
tures in the lens galaxy halo. 

However, these effects will strongly depend on the sub- 
halo mass and density profiles. In paper I, we have shown 
that the image separation caused by substructures in the 

' sec CASTLES website: |http: //cfa- www.harvard.edu/castles/ 
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Fig. 7. — Average number of substructures projected on an 
extended source with radius r s = 1 kpc at z s = 2 as a function 
of projected radius. For a host halo at z; = 0.5 with mass M = 
1.8 X 10 12 M Q we vary the minimum subhalo masses with m m j n 
. x 1O 6 M (dotted line), 1O 8 M (dashed line) and 1O 9 M (solid 
line), respectively. We assume a maximum magnification factor of 
30 at the Einstein radius of the host halo. 

mass range < 10 10 M Q will only be resolvable if one can 
assume a steep inner density profile similar to SIS and 
M99. Even for these favorable density profiles, there will 
be a minimum mass for which the image separation drops 
below the resolution of any present or planned observa- 
tional facilities. For instance, for the currently available 
EVN 8 , the minimum subhalo mass which could be re- 
solved under optimal conditions is approximately 4 x 10 6 
M© in case of an SIS profile but increases to 3 x 10 7 M© 
for an M99 profile. Figure [7| shows how the expected 
number of intervening substructures N depends on the 
minimum subhalo mass m m i n . For a host halo of mass 
M = 1.8 x 10 12 Mq at zi = 0.5 and a source of size r s = 1 
kpc at z s = 2, we compute TV as a function of projected 
radius from the host halo center for minimum subhalo 
mass m min = 4 x 10 6 M Q , 10 8 M Q and 10 9 M Q , respec- 
tively. Since the subhalo mass function predicts about 
equal mass within each logarithmic mass bin, most sub- 
halos will be of low mass (n(m) oc rn~ ' 9 ). Thus, the ex- 
pected number of substructures projected on the source 

8 www.evlbi.org 



is very sensitive to the minimum subhalo mass that can 
be resolved. 

Therefore, angular resolution will be crucial when at- 
tempting to detect CDM substructure via its lensing 
effects on background quasars and submilliarcsecond- 
resolution facilities will be required. We assess that 
the upcoming VSOP-2 satellite may provide the best 
prospects for such a detection. However, one must be 
careful not to confuse internal structures found in quasars 
observed at radio wavelengths with subhalo lensing sig- 
nals. Quasars already macrolensed on arcsecond scales 
can be used to test that it is possible to distinguish 
between the two, since internal structures should be 
mapped in all of the macro-images while subhalo lens- 
ing will only affect one of the images. Therefore, such 
system will be ideal to train this technique of identify- 
ing the strong lensing signal from dark matter substruc- 
tures. However, even single-imaged systems with rather 
large projected lens-source distance could be used to set 
upper limits on the abundance and nature of subhalos 
depending on the amount of distortion signals detected. 
Taking the above into account, one may be able to use 
this technique to put constraints on dark matter subhalos 
predicted by simulations in the near future. 
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